KINETIC MODELS FOR POLYMERS WITH INERTIAL EFFECTS 



PIERRE DEGOND AND HAILIANG LIU 



Abstract. Novel kinetic models for both Dumbbell-like and rigid-rod like poly- 
mers are derived, based on the probability distribution function /(t, x,n,h) for 
a polymer molecule positioned at x to be oriented along direction n while em- 
bedded in a h environment created by inertial effects. It is shown that the 
probability distribution function of the extended model, when converging, will 
lead to well accepted kinetic models when inertial effects arc ignored such as 
the Doi models for rod like polymers, and the Finitely Extensible Non-linear 
Elastic (FENE) models for Dumbbell like polymers. 
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1. Introduction 

In this paper we derive novel kinetic models for both Dumbbell like and rigid-rod 
like polymers in the presence of inertial forces. The model is to describe dynamics 
of the probability distribution function embedded in high dimensional configuration 
space due to inertial effects. Wc then prove that the limit equation of the new model 
when inertial force vanishes leads to current models with no inertial effects such as 
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the FENE model and the Doi model, respectively. This illustrates consistency of 
our kinetic models with existing models. 

The wide range of applications of polymer materials has attracted new areas of 
academic and industrial research. The synthesis of different type of polymers has 
enlarged the range of applications of polymer materials to areas where mechanical 
properties are important. Materials made up of macro-molecules such as polymers 
display properties that completely differ from those made from small molecules. The 
description of polymer dynamics is often based on large assemblies of molecules, the 
characteristics could be modeled in terms of their statistical properties. 

Most polymers are long chains or branches of repeated chemical units. The full 
description of each atom in the polymer by molecular dynamics is not feasible for the 
huge computational effort. Coarse grained models are often expected with macro- 
scopic space and time properties of complex fluids. Typical models such as bead- 
spring chain for flexible polymers and the rigid rod model for liquid crystalline 
polymers have been established by the pioneers in polymer science. 

In general, the flow modeling of polymers has to take into account the internal 
structure, characterized by both positional and oricntational order of phases. Such 
incorporation is often done by adding new balance equations to those that govern 
structure-less Newtonian fluids. These new balances must be evaluated from the 
behavior of polymers. According to the relative size of the bending persistence size 
and the length of the polymer, two canonical types of polymers are widely studied: 
the Dumbbell model and the rigid-rod like model. In modeling motion of polymers, 
it is essential to explore an accurate method of solutions of the Langevin equation 
for particles undergoing Brownian movement (rotational or translational) under the 
influence of external fields. A large number of reviews, text books and monographs 
on the theory, applications and rhcology of polymeric materials have appeared in 
the literature, see e.g. [IHl HS1 12DI H3 HDl HH ESI HI ESI 1231 EH] - 

There are three main levels of description of polymeric fluids: atomistic modeling, 
kinetic modeling [III 120] , and the macroscopic approach of continuum mechanics 
[36] . We shall exploit the kinetic approach. Models of kinetic theory provide a 
coarse-grained description of molecular configurations wherein atomistic processes 
are ignored altogether (Doi and Edwards [TTj . Bird et al [36], and Ottinger[33]). 
Kinetic theory models for polymer solutions are most naturally exploited numer- 
ically by means of stochastic simulation or Brownian dynamics methods [33]. A 
kinetic theory model when equipped with an expression relating stress to molec- 
ular configurations plays an important role in developing micro-macro methods of 
computational rhcology [35[ I23j. In current kinetic theory models for polymers, the 
inertia of molecules is often neglected. However, neglect of inertia in some cases leads 
to incorrect predictions of the behavior of polymers. The forgoing considerations in- 
dicate that the inertial effects are of importance in practical applications, e.g., for 
short time characteristics of materials based on the relevant underlying phenomena. 

It is thus the goal of this paper to model dynamics of the density distribution of 
polymers when the inertial force is no longer ignorable. More precisely we shall be 
particularly interested in modeling two canonical types of polymers: Dumbbell-like 
and rod-like polymers, which when inertial forces are not considered have been well 
understood. We first derive kinetic models including inertial effects from particle 
dynamics (continuum limit in the Brownian motion), we then show the limit of the 
augmented models when inertial forces vanish leads to the inertia-free model. 

We now summarize our main results for two types of polymers. 



3 



1.1. Dumbbell-like polymers. A mcromolcculc is idealized as an 'elastic dumb- 
bell' consisting of two 'beads' joined by a spring which can be modeled by an end- 
to-end vector n. Here n is in a bounded ball B(0, no), which means that the exten- 
sibility of the polymers is finite. Let f(t,x,n,p,q) denote the distribution function 
of Dumbbell-like polymers on the space variables x G K d , d = 2, 3, the translation^ 
velocity p G M. d , the end-to-end vector n as well as the oricntational velocity q G M d . 
And t is the time. The novel kinetic model to be derived is 

dtf + V x ■ (pf) + V„ ■ (qf) + V p ■ (-^{p- u(x))f) (1.1) 

m 

+ V q -((--{q-n-VM*))-—) f) =^p[A p / + A g /], 
\\ m m J J m z 

where £ is the frictional coefficient for the beads with mass m, k B is the usual 
Boltzmann constant, T is the absolute temperature, and F is the spring force be- 
tween beads. The force usually derives from a potential, and has different forms for 
different models. For the well known FENE potential 

F Hn 



where H is a spring constant |36j . 
The above model under the scaling 

™ _ ,2 



'mp~^p, \jmq—tq, 
leads to 

e 2 d t f + eV x ■ (pf) + eV„ ■ {qf) 

+ eV p • (£u(x)f) + eV, • [(Cn • V*u - 2F)f] = Q{f), (1.2) 

where 

Q(f) ■= CV P • (pf + 2k B TV p f) + CV 9 • (qf + 2k B TV q f). 
Our result for Dumbbell like polymers could thus read as follows: 

Theorem 1.1. The limit e — > of f is given by f° = pM where p = p(t, x,n) > 
and M are given by 



p(t,x,n) = / f°(t,x,n,p,q)dpdq, M = exp 



p 2 + q 2 
' ^k B T 

Furthermore, p(t, x, n) satisfies the following kinetic equation 

2i*\ 2k B T 2k B T 
o t p + V x ■ (u(x)p) + V„ ■ ((n • V x u —)p) = — — A x p H —A n p. (1.3) 

1.2. Rod-like polymers. Polymers are idealized as rods of fixed length. The ori- 
entation space is n G S d ^ 1 . Let f(t,x,n,p,uj) denote the distribution function of 
rod- like polymers on the space variables x G R d , the translational velocity pel 11 , 
the orientational vector n as well as the angular velocity uj. Here lu is on the tangent 
bundle T n S . And t is the time. Our rescaled kinetic model can be formulated as 

ed t f + V, • (pf) + TZ ■ (uf) + C*V P • (u(x)f) 

+ V w • ((( r n xV x u-n-K- U)f) = -Q(f), (1.4) 

e 

where 

Q(f) = CtV p • (pf + k B TV p f) + C r V w • (fuj + fe B TV w /). 
Here e denotes the inertial parameter similar to (|1.2|) , u is the fluid velocity, and U is 
certain interaction potential of rods. 1Z — n x V„ is the rotational gradient operator, 
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and ks,T denote the Boltzmann constant and absolute temperature, respectively. 
Ct, £ r are frictional coefficients in x and n directions. 
Our result for rod-like polymers then reads as follows: 

Theorem 1.2. The formal limit e — ► of f is given by f = pM where p = 
p(t, x, n) > and M are given by 

p(t,x,n) = J f°(t,x,n,p,q)dpd n uj, M = exp ^ ■ 

Furthermore, p(t, x, n) satisfies the kinetic equation 



dtp + V K • (u(x)p) +TZ-(nx W x u ■ p) = D t A x p + D r TZ ■ 



Tip 



1IU 



(1.5) 



where 



k B T 

Ct ' 



D r = 



k B T 
Cr ' 



Our derivation of these kinetic models is based on establishing motion laws of 
polymer molecules, followed by a conversion into the kinetic description. The formal 
limit when inertia vanishes is justified by taking the classical approach for hydro- 
dynamic limits. To this end a rcscalling is adopted, so that the collision operator 
is set on fast scale, and the dissipation of the collision operator drives the states to 
the unique equilibrium states M. Derivation of models for Dumbbell- like polymers 
and the formal limit justification are given in §3, and those for rod-like polymers are 
given in §4. Some concluding remarks are presented in §5. 

Finally, we wish to close this section by pointing to a vast body of recent work 
on mathematical treatment of kinetic theory models for polymers, their constitutive 
models as well as their coupling with fluid models (so called micro- macro models), 
see e.g. [3Zll[Illl[3[I3[22[Il[l[Zl[^ 

l32l [34l |4T1 [40] and references therein. These should provide guidance in establishing 
the various levels of mathematical theory of kinetic models developed in this work. 



2. Kinetic models for Dumbbell polymers 

2.1. Equations of motion with non-trivial inertia. We consider a polymer 
consisting of two-beads connected by one-spring. Each bead as a coarse grained 
particle represents several chemical units and experiences four kinds of forces in the 
dilute case where there is no interaction for intcr-and intra-dumbbclls. 

Assume the two beads are positioned at x\ and x 2 , the Langevin equation of the 
beads are just balance of different forces expressed as 

m\±\ = -Ci(ii - u{x\)) + Fx + axWx, 

m 2 x 2 = -(2(x2 - u(x 2 )) + F 2 + a 2 W 2 . 

Here Wi are independent standard Brownian motions, C,i are frictional coefficients 
for the i-th bead with mass mi, Oi = ^/SfcgT^ by fluctuation-dissipation theorem, 
where ks is the Boltzmann constant, and T is the temperature. The spring forces 
Fx = —F 2 = F by Newton's Third Law depends only on the end-to-end vector n. 

A natural configuration for this underlying polymer includes both position of the 
polymer with these two beads connected to one spring 

Xx + x 2 



and the end-to-end vector denoted by n — x 2 — xx- We assume that mi = m 2 = m 
and (i = ( 2 = C- Then the difference and average of the above two equations gives 
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a coupled system 



mi = -C{x - u{x)) + y/4k B T(Wi, 



) - 2F + V4^TCW2 



(2.1) 
(2.2) 



The approximation utilizes the equation u(x) ~ (w(xi) + u(x2))/2, u(x2) — u(xi) ~ 
V^u • n, and the two new Brownian motions (W2 ± Wi)/V% (still denoted by 

If the inertia were ignored for this two beads spring system, we would have the 
following stochastic equations (SDE), 



u(x) 



l4k B T 



Wi, 



2F Uk B T ■ 
n = n-V x u — + W — —W 2 . 



(2.3) 
(2.4) 



The differential operator or the infinitesimal generator corresponding to the process 
(x(t),n(t)) is given by 



2F 2knT 
L = (u{x),V x ) + (n ■ \7 x u - —, V n > + — —(&x + An), 



(2.5) 



where 



d_ 

d Xl 



' d„ 



with (, ) denoting the Euclidean inner product in R d . The configuration of Dumbbell 
like polymers can be illustrated in Figure [T] 




r-*» n, 



Figure 1 . A Dumbbell like polymer 



^Here we use the fact that for any a G [0,1], aWi(i) ± s/l — a 2 W2(t) remains a Brownian 
motion. 
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When incrtial force is significant, one has to take the inertial effects into consid- 
eration. We now rewrite the system (|2.f l) - (|2.2[) into a first-order system 

x=p, 
h = q, 

p = (p - u(x)) H VKi, 

m fa 



9 = — (q- n - y x u(x)) 1 Pv 2 - 

m mm 
This system can be considered as a degenerate stochastic differential equation with 
a singular diffusion matrix of Ad x Ad. The differential operator corresponding to the 
process (x(t) , n(i) , p(t) , q(i)) is also uniquely defined. In fact for the above dissipa- 
tive stochastic differential equation with well-defined initial data, a unique solution 
(x, n,p, q)(t) up to an explosion time is ensured if u(x) and F are smooth functions 
of the configuration variables. 

2.2. Kinetic description. We now formulate the kinetic description of the above 
motion laws. For a given instant t, the random variables X = (x,n,p,q) can be 
characterized by a probability density function (PDF) f(t,x,n;p,q) defined by 

f(t, x, n; p, q, )dX = P[{X < X < X + dX}] . 

In words, the right-hand side is the probability that the random variable X falls 
between the sample space values X and X + dX for different realizations of the 
polymer motion. Thus the above statistical ODE system can be converted into a 
PDE of the form 

dtf + V, • (pf) + V„ ■ (qf) + V p ■ (-^(p - u(x))fj (2.6) 

+ V 9 ■ ((--(q - n ■ V x «(a:)) -—)f)= ^^[A p f + A g f). 
\\ m m J J m 

It is known that in the case with no inertial forces the kinetic equation is also of the 
Fokker-Planck type. The model follows from the motion law (|2.3[) - (|2.4[) as follows: 

[n-V x u- — \p\ =— ^—A x p+— ^~&nP, (2.7) 

where p(t, x, n) denotes the corresponding probability density function. 

Remark 2.1. The polymer scale is much smaller than the fluid scale, the diffusion 
in space is often ignorable. Here we make no such distinction. 

It is natural to ask whether the effective dynamics in (|2.6[) is dominated by (|2.7p 
when inertia vanishes. As a first step we shall identify an effective equation satisfied 
by the limit of / when m | 0, if the limit exists. In the remainder, we will always 
suppose that / is as smooth as needed. 

2.3. Scaling. We make the following scaling 

m = e 2 , 

to obtain 



mp — > p, \/mq — > q, 



e 2 d t f + eV x ■ (pf) + eV„ ■ (qf) + eV p • (Cu(x)f) + eV q ■ {((n ■ V x u - 2F)f] = Q(f), 

(2.8) 

where 

Qtf) := CV P ■ (pf + 2k B TV p f) + CV, • (qf + 2k B TV q f). 
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The remaining task in this section is to justify the following 

Theorem 2.1. The limit e j of the / £ is given by f° = p(t, x, n)M(p, q), where 

satisfying Q(M) = and p(t,x,n) solves the kinetic equation {2.7]) . 

As usual the limit equation that is associated with (|2.8p does not depend on 
details of the operator Q. Rather, it depends on Q possessing certain properties 
related to conservation, dissipation, and equilibria that are stated below. Note the 
operator Q acts on the augmented variables (p, q) only and leaves other variables 
(t, x, n) as parameters. Thus we list some properties of Q as an operator acting on 
functions of (p, q) only, which are essential in the derivation of the limit equation. 



2.4. Properties of Q. Note that the operator Q(f) has the following properties 
• The operator can be written as a conservative form 



Q(f) = 2(k B TV p , q • (MV^ 
The conservation form leads to 

Q(f)dpdq = 



for every / € L l (dpdq). This relation expresses the physical laws of mass 
conservation during the action of Q. 

Moreover, this is the only such conservation law. This means that 



J Q{f)ipdpdq = 0, V/ & tp G span{l}. 



Dissipation. There is a nonnegative function rj(f) that is an entropy for the 
operator Q. This means that for rj(f) G ^(dpdq) we have 

v'ifMf) < o, 

whose vanishing characterizes the local equilibrium of Q. For flexible poly- 
mers, Tj'(f) = f/M since 

D Q {f) - J Q(f)jjdpdq = - J M (v p , q (j^yj dpdq < 0. 

• Existence of equilibria. From the above dissipation property we see that the 
function /(p, q) such that Q(f) = forms a space identified through 

/ 

— G span{l}. 

This equilibrium is unique up-to a constant multiplier. We normalize the 
density function then 



feq — 



Mdpdq 



M. 
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2.5. Asymptotic limit as e J, 0. Equipped with these properties we now investigate 
the asymptotic limit. Define 



p e (t,x) = J fdpdq, 
J{{t,x) = i J pfdpdq, 
J£(t,x) = ^ /?. 



fdpdq. 

Then integration of (|2.8|) together with conservation property of Q yields 

d t p e + V x - Ji+V„- J 2 = 0, Ve>0. (2.10) 

We suppose that as e | the limit of / is identified by f°. It follows from the 
equation (|2.8p that 

f°&Q-\0). 

This implies that /° = M(p, q)p°(t, x, n), for the operator Q does not act on (t, x, n). 
Hence the probability density function / takes the following form 

f e (t,x,n;p,q) = p t (t,x,n)M(p,q) + eg e , 

where 

\imp e (t, x, n) = p (t, x, n). 
Note that / pMdpdq = J qMdpdq = 0, we thus have 

J{ = \j Pfdpdq = J pg'dpdq, (2.11) 

J 2 = \j ifdpdq = I qg'dpdq. (2.12) 

The limit of these flux functions depends only on the limit of g e . 
Linearity of the operator Q and Q(M) = leads to 

Q{f c ) - Q(p'(t,x)M(p, q) + eg') = eQ( 5 e ). 

Using the equation (|2.8|) we obtain 

Q(g e ) = -Qin 

= ed t r + V x ■ (pD + V„ ■ (qf e ) + V p • ((u(x)r) + V q ■ ((Cn • V x u - 2F)n 
= V x ■ (pMp e ) + V„ • (qMp e ) + V p ■ ((u(x)p e M) 
+ V q ■ ((Cn ■ V x u - 2F)Mp e ) + 0(e). 

Let g e — > g° as e — > 0, then it is clear that 

Q(g°) = P ■ VxP°M + q ■ V n p°M + Cp°u(x) ■ V P M + p°((n ■ V x u - 2F) ■ V q M. 

Again note that Q is a linear operator, we may assume the ansatz for g° as follows. 

.9° = a ■ V x p° + b ■ V n( o° + (p°c ■ u(x) + p a ((n ■ V x u - 2F) ■ d. 

In comparison with the expression of Q(g°) above, we see that it is sufficient to find 
a, b, c, d such as 

Q(a) = pAI, 
Q(b) = qM, 
Q(c) = V P M, 
Q(d) = V q M. 







These relations arc understood in the sense that the operator acts on each component 
of the underlying vector. In order to identity (a, b, c, d) we state the following lemma: 

Lemma 2.1. Let Q be the collision operator and g G L 2 (R P x R q ) such that 
J gdpdq = 0. The problem 

QW=g, (2-13) 

has a unique weak solution in A = {tp, J ipdpdq = 0, tp € H 1 }. 

This lemma, following from applying the Lax-Milgram theorem to the corre- 
sponding variational formulation, ensures that there exists unique a, b, c, d in the 
space {0| / cj)dpdq = 0} since 



J {pM, qM, V P M, \7 q M}dpdq = 0. 



We also note that from 
it follows 

a b 
d = — - 



2k B T 7 2k B T' 
It remains to determine only a and b. Let a = aMp we have 

Q(a) = -a(pM = pM, 

which gives a — — ^, leading to a = —^pM. Similarly we obtain b = —^qM. 
Clearly we see that 

J {b, d} ■ pdpdq = 0, J {a, c} • qdpdq = 0. 

These enable us to determine the asymptotic limits of the fluxes. 

J° = J pg°dpdq 

(a ■ V xP °)p + f(b- V n p°)p + Cp° I (c • u(x))p + p° f d ■ (Cn ■ V x u - 2F)p 



V^-^-Y-^lMdpdq. 



2k B T ) C 

4 = J qg°dpdq 

= J (a ■ V x p°)q + jb- V nP °q + Cp° J (c • u{x))q + p° J d ■ ((n ■ V x u - 2F)q 



= Jb- (v nP ° --^(Cn- V x u - 2F)^j qdpdq 

V„p° - ■ V x u - 2F)\ J Mdpdq. 



Now the equation (|2.10p divided by J Mdpdq asymptotically converges to 
d tP ° + V x ■ (u(x)p°) + V„ • ((n • V x u - ^V) = ^A xP ° + ^ 
This is exactly the equation (12. 7p derived from ignoring inertial forces. 
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Remark 2.2. With these kinetic equations for Dumbell like polymers, it is natural 
to understand how these polymers contribute to the macroscopic flow governed by 
a coupled Navier-Stokes system 

d t u + (u ■ \7 x u)u + \7 xP = —A x u + X7 X ■ r, V x ■ u = 0, 

Re ReUe 

with the stress determined by 

r = J (n (g) F) f dndpdq. 

Here -De is called Deborah number, which is the most important parameter in non- 
Newtonian fluids. Re and 7 are the Reynolds number and viscosity ratio, respec- 
tively. The tensor force follows the case when inertial force is ignored, the derivation 
of the tensor force with inclusion of inertial effects is beyond the scope of this work. 

Remark 2.3. The spring force is determined by a potential function U through 
F = —V n U. Different potential leads to different models. Two choices are com- 
monly used: 





Force 


Potential 


Hookean spring 


Hn 


\Hn 2 


FENE spring 


Hn 


-\Hnlln(l-{n/n f) 


l—(n/n ) 2 



where no is the maximum extension of the beads connector. 



3. Rigid rod-like polymers 

Though many polymers are flexible, there is still a large class of polymers which 
are not flexible and assume a rod-like structure. Rod-like polymers have some pe- 
culiar properties and have attracted a great deal of attention. 

We consider rod-like molecules in concentrated regime. Rod-like polymers can 
have only two kinds of motion, i.e., translation and rotation. The translational 
Brownian motion is the random motion of the position vector x of the center of 
mass, and the rotational Brownian motion is the random motion of the unit vector 
n (\n\ = 1) which is parallel to the polymer. We shall build a kinetic model for the 
probability distribution of oricntational motions of rod in every point of phase space 
(x,p). This serves as a microscopic equation, which is expected to be coupled with 
the macroscopic equation (the Navier-Stokes equation) for the fluid velocity. 

For the convenience of calculations in what follows, we introduce a local coordinate 
on the sphere § 2 as q — (6,<fi), and set n = (sinBcoscj), sinOsiruf), cos6) with 6 € 
[0, 7r) and <fi € [0, 2tt). The unit vector eg = (cos0cx)S(f>,cos6sin(f>,—sin6) and = 
(— sincf), cos(f>,0) . We note that ^ = eg, ^ = sinQe^. Thus, any tangent vector on 
the sphere is written as 

n = 9eg + (jisinOe^, 

which gives 

u) = n x n = Qe^ — sin9<j)eg. 

The gradient is 

V n = egdg + -^tA- 

The rotational gradient is 

cos(j)cos0 . cos6sin(j) \ T 
sinO — * ~ sm ^" e ' ~n — "<t> + c.os<pdg,d^ I . (3.1) 



1Z = n x V„ = 
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Divergence of a vector A = A g e g + A^e^ is 

V„ ■ A = -^—de{sinQAg) + -^deA^. 
sina stnt) 

The gradient in terms of lj is defined only for a given n, and will be understood from 
now on as 

n ■ V u f = 0, V/. 

To be more specific, we regard the identical liquid crystal molecules as inflexible 
rods of a thickness b which is much smaller than their length L, as illustrated in 
Figured 



u 




b 



Figure 2. A rod like polymer 



3.1. Translational Brownian motion. Let U be an interaction potential, the 
force thus induced is —\7 X U. The motion law is 

mx = -( t (x - u(x)) - V X U, 

where is the friction coefficient and u(x) is the fluid velocity field. If / is the prob- 
ability distribution in x, the Brownian force is expressed as Fb = —Vx{kBTlnf), 
where —ksTlnf is the chemical potential. For nontrivial incrtial force, the distri- 
bution needs to be accounted in an extended environment with inclusion of p — x. 
The corresponding translational Brownian force thus reduces to 

F B =V P W, W = -^^lnf. (3.2) 
m 

This can be justified by a similar derivation based on Brownian motions as that in 
Section 2. The scaled coefficient reflects balances between the friction force and the 
inertial force. Putting together we have the following translational motion law 

x=p, (3.3) 

V = --0 - u{x)) - -V X U - ^-V p (k B Tlnf). (3.4) 
m m m z 

In what follows we shall conveniently use the chemical potential to describe the 
Brownian force. 
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3.2. Rotational Brownian motion. We consider a rod rotating with angular ve- 
locity uj, then the rotational motion can be described as 

[J]uj = T, 

where J is the moment of inertia and T is the total torque. 

1. Rotational frictional force. 

Consider a rod of length L placed in a viscous fluid with fluid velocity u(x), where 
x is the center of mass of the rod. 

The rod is parameterized by s, ranging from — L/2 to L/2, then the position 
vector of the s-point on this rod is written as 

x(s)=x + ns 7 x(0)=x, -L/2 < s < L/2. 

Let v(s) and F(s) be the velocity of this point, and forces acting on it. The velocity 
v(s) is expressed by the angular velocity uj 

v(s) = uj x (x(s) — x) = uj x (ns). 

The frictional force at x(s) is 

F(s) = -Z(s)(v(s)-u(x(s))), 

where £(s) is the frictional coefficient, being symmetric £(s) = £,(~s). Note that 

u(x(s)) ~ u(x) + \7 x u(x(s) — x(0)) = u(x) + sV x u • n. 

The frictional force thus reduces to 

F(s) = —s£(s)(lo x n — V x u ■ n) + £(s)u(x). 

Thus the total torque induced by the frictional force acting on the rod is 

i r L ' 2 

T = - (x(s) - x) x F(s)ds 

L J-L/2 

— —£,2n x x n — V x w • n) + x u(x), 

where 

1 f L/2 

tk = T s k ti{s)ds, k = l,2. 

L J-L/2 

Symmetry of £(s) = ^(— s) leads to = 0. Then using n ■ u> = we have 

T = — C, r [(Id — n ® n)uj — n x (V x u ■ n)] = —£ r [uj — n x (V x u ■ n)] , 

where £r = & denotes the rotational friction coefficient. 

2. Thermodynamic potential force. 

From [J]uj =T it follows 

±E = T. U , E=y { J)uj. 

Let the potential be denoted by U. Then d(E+U) = gives dE = —dU = —V n U-dn. 
Note that from |n| = 1 we have 

uj = n x n, 

which leads to 

uidt = n x dn. 

We thus have 

T-(nxdn)=T-ujdt = dE= -V n U ■ dn. 
Therefore T X n= — V n J7, i.e., 

T = -KU, (3.5) 
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where 1Z is the rotational gradient operator given by (|3.1[) . 

3. Rotational Brownian force. 
The following property will be used on our derivation of forces below. 

Lemma 3.1. For a fixed vector a, let x — a x y. Then for any function smooth g, 

Vy.g = -ax V x g 

and 

a ■ Vj,<? = 0, Vg. 

Proof. Let as the usual permutation symbol, then Xi = ^ . k CijkOjyk- We thus 
have 

(V v flOm = ^ d x% g^p- = ^<V r //^ /f" 1 '/-./. 

Oym ■ ■ , 

I " % J,k 

ij hi 

This gives the desired relation. We now show the second claim. 

^ Q"m&y m g ^ ^mijdxi9CL 7n Oj ^ ^ fmij^i5%^mi 

which ensures that a ■ d y f = holds for any smooth function g 

□ 

Consider a rod again with center of mass at x, and momentum p = at rest, the 
point vector of rod is x(s) = x + sn with s ranging in [—L/2, L/2]. Correspondingly 

p(s) = sh = suj x n. 

From this we obtain that 

oj = —n x p(s). (3-6) 
s 

Also rod symmetry implies that 

p(-s) = -p(s). (3.7) 
For any fixed s the associated Brownian force is calculated by 

F B (s) = V p(s) W, W = -^^-lnf. (3.8) 

m 

Note that for each fixed s and vector n, the mapping from p(s) to lo is well-defined. 
The result in Lemma 13-11 gives 

F B {s) = V p{s) W = --nx\7 u W. 
Thus the corresponding Brownian torque is determined by 

-L/2 



i r /z 

T B = y I snx F B(s)ds, 

L J-L/2 
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which when combined with the above calculations leads to 

1 f h l 2 ( 1 \ 

T B = — I snx — n x V^W ds 

LJ-L/2 \ S J 

= -n x (n x V W W) 
= (Id-n® n)\7 u W 

= -^^-{Id - n ® nWJlnf) 
m 

( r k B T 



-V u (lnf). 



in 

Note that the moment of inertia for the rod is 

mL 2 



[J]= j(Id- n®n), j = 



12 

Using Co ■ n = 0, we thus have 

[J] ■ d) = jd). 

This together with all forces involved leads to the following motion law 

h = oj X n, (3.9) 
^-Co = -C r (w - n x (V B u) • n) - ^^-V^lnf) - TZU, (3.10) 

12 TO 

where ?7 is the interaction potential. 

Remark 3.1. The potential can be defined as a generalized Onsager's potential: 

U = v 2 k B TbL 2 / / / B(x,p;x',p')\nxn'\f(t,x',p',n',uj)dn'dx'dp'd n uj. 

J JR P J\n'\ = l Jn-a>=0 

Here i? x and R p is the configuration space for variables x and p. B is some localized 
kernel. 

3.3. Kinetic equations. We shall derive a kinetic equation in the phase space 
(x,p,n,uj) with (x,p) e R 3 x I 3 and (n, w) € T§ 2 , where 

TS 2 := {(n,u>)\ neS 2 , lu g T„§ 2 }, (3.11) 

which is usually called the tangent bundle of the manifold § 2 . We start from the 
continuity equation of a formal form 

dtf + y x ■ (if) + V P ■ (pf) + V„ • (nf) + V A • (nf) - 0. (3.12) 

This equation can be simplified when restricted on the tangent bundle: (n, u>) £ TE> 2 . 
First we state the following 

Lemma 3.2. Let (n,oS) be any element of the tangent bundle TS 2 , mapped to an- 
other element (m, n) <G T§ 2 such that 

to = n, h = lo x n, 

then for any smooth function f we have 

V n / = V m / - uj x V A /, V w / = n x V*/. (3.13) 
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Proof. Let tijf, be the usual permutation symbol, then nj = J2jk e ijkUjn k . Further 
we have 

i j,k 

where (5^ is the Kronecker delta symbol. This gives the first relation in (|3.13[) . The 
second relation follows from 

i j.k 



Lk 



□ 



Using similar arguments we have the following: 
Lemma 3.3. Let g be any smooth vector function, then 

(n x V„) -g = -V„ ■ (n x g), (n x V w ) • g = -Vu, ■ (n x p). (3.14) 
Note from uj = n x h we have w = n x fi leading to w ■ n = 0. Also 

n = u> x n, h = Co x n — |w| 2 n. 

Equipped with above lemmas and relations, we are able to reduce the equation 
(|3.12p . We only check the last two terms on the left-hand of (|3.12[) . First, 

V„ • (nf) = V„ -(uxnf)=n- (vf), 

where 1Z := n x V„ is the rotational gradient operator. For any smooth function g, 
Lemma 3.2 and Lemma 3.3 imply that 

= — (n x W u )g. 
This enables us to simplify the last term in equation p. 121) 

V h ■ (nf) = -{n x V w ) -{{uxn - \uj\ 2 n)f) 
= V w ■ (n x (w x n — |w| 2 n)/) 
= V w • ({Id- n®n)iof) 
= V w ■ (cbf). 
Therefore the effective kinetic equation becomes 

dtf + V, • (p/) + V p • (p/) + ft • (a;/) + V w • (u//) = 0, (3.15) 

where 

P = — — [p - u(x)) V X U - -^rV p (k B Tlnf), 

to to m z 

Co = — — (u! — n x V x u(x) ■ n) 1ZU — ^rV w (/csTZn/). 

TO TO TO^ 

Note here the coefficient L 2 /12 has been absorbed in both the £ r and the potential 
U, which is independent of both p and to. 
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3.4. Scaling. We are interested in the solution behavior when inertia vanishes. 
We now make the following scaling 

to = e 2 , tp — > p, euj^uj, 
under this scaling the system (|3 . 1 5f> written into the new variables becomes 
ed t f + V x ■ (pf) + K ■ (w/) + ( t V p ■ (u(x)f) 

+ V w • ((C r n x V x u ■ n - K ■ U)f) = \q(J), (3.16) 

where 

Q(f) = C t v p ■ (pf + k B rv p f) + c,v w • {fu + fc B rv w /). 

Our next task is to investigate the formal limit e — > of this problem. 
Theorem 3.1. The limit e J, of the / £ is given by f = p(t, x, ri)M(p,co), where 

"-"(-isf) w 

satisfying Q(M) = and p(t,x,n) solves the following kinetic eguation 



KB-l 



(3.18) 



dtp + V K • (u(x)p) +K-(nx V x u ■ p) = D t A x p + D r K 
where 

_ k B T _ k B T 

Remark 3.2. (i) When the coefficient D t = 0, the equation (|3.18[) is the Doi's kinetic 
model for rod-like polymers with inertial forces ignored. This equation is also called 
the Smouluchowski equation in literature. 

(ii) If the translation diffusion has different strengths and U also depends on x, one 
needs to change DtA x p to 

V x ■ \ [D ll n®n + D ± (Id-n®n)][V x p+ T ? 7 =V x U] 
L k B l 

The above justification procedure remains valid. 

3.5. Properties of Q. The operator Q here acts on the augmented variables (p, ui) 
only and leaves other variables (t, x, n) as parameters. Thus we list some properties 
of Q as an operator acting on functions of (p, ui) only, which are essential in the 
derivation of the limit equation. 
Set 

/ p 2 + U? * 

with uj satisfying uj ■ n = for any fixed n £ § 2 . 
This form enables us to conclude the following 

Lemma 3.4. The operator Q has the following properties: 
(i) The operator Q can be written as 
Q(f) := Ctk B TV p ■ \ MV P f-^J J + ( r k B TV u ■ [MW^ (Jj)) ■ 
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(ii) Conservation. For every f G L 1 (dpd n uj) , 

Q(f)dpd n uj = 0. 



Moreover, this is the only such conservation laws. In other words, for any 
f G ^(dpdnuj), 

Q(f)ipdpd n u = 0, V/ V> e span{l}. 

(iii) Dissipation. There is a nonnegative junction rj(f) that is an entropy for the 
operator Q. This means that for r)(f) G L l {dpd n ui) we have 



D Q (f) = J Q(fW(f)d P d n oj<0. 
In the case of rod-like polymers, rj'(f) = f/M and the dissipation production 



is 



-k H T 



C t J M (v p dpd n iu + Cr J MV W (J^j dpd n L 



< 0. 



(iv) Existence of equilibria. From the above dissipation property we see that the 
function f(p,u>) such that Q(f) = if and only if 

f 

— G span{l}. 

Proof of Lemma 13.41 (i) The conservation follows from a direct calculation, (ii) 
The conservation form in (i) leads to the null integration. Since the J Q(f)ipdpd n oj = 
holds true for any /, we let / = ipM to have 



Q(ipM)ipdpd n u> = 0. 
The left-hand side of this relation is non-positive as 



= -k R T 



Ct / M (V P V) dpd n uj + Cr M{Wu,ip) 2 dpdnUJ 



< 0. 



Thus i/j must satisfy both V p -0 = and W^ip = 0. Then one has n ■ = 0. These 
ensure that ip must be independent of (p, oS). 

(iii) Dissipation property comes from integration by parts on the tangent plane. 

(iv) From the dissipation inequality we see that if Q(f) = 0, then D®(f) = 0. 
The non-negativity of yields 

/ 

— G spanjl}. 

3.6. Asymptotic limit as e j 0. We now investigate the formal limit of the prob- 
lem, assuming all involved functions are smooth and convergence holds true as 
needed. 

We suppose that / — > f° as e — > 0. Then, from the scaled kinetic equation (13.161) . 
Q[f] = O(e) and we deduce that Q(/°) = 0. 

By property (ii) we have f° = pM, with p = p(t,x,n) > and n G § 2 . Note 
that Q acts only on (jp, uS), p arc functions of (t,x,n). To find this dependence, we 
use the generalized collisional invariants. We integrate the equation with respect to 
(p,u>) to find the continuity equation 

d t p e + V x ■ J{ + K ■ Jt, = 0, Ve > 0, 
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where the density and fluxes are defined by 

p e (t,x,n) = J fdpd n uj, 

J{{t,x,ri) = i J pfdpd n u; 7 

JJ(i, X, n) — — J uufdpd n uj. 

Here the integration over w is interpreted as for any fixed n with (n, cj) G T8 2 . In 
the limit e -> 0, p £ -> Cp° and Jf -> CJ° (i = 1, 2) with C = / Mdpd n u, and we 
obtain 

<9tP° + V x • J? + TZ ■ J 2 ° = 0. (3.20) 

Thus we may assume that 

f = p e (t,x,n)M(p,u) + eg e (t,x,p,n,uj), 

which leads to 

Q(f e ) = zQ(g e )- 

To determine the limiting flux we need to explore the limit of g e . Using the kinetic 
equation (|3.16[) we have 

Q(g e ) = -Q{f) 

= ed t f + V x ■ (pf) + K ■ (luT) + V p • {Ctu{x)n + V w • ((C r n • V x u - KU)n 
= V x ■ (pMp e ) + K ■ {uMp e ) + V p ■ ((tu(x)fM) 
+ V w • ((C r n • V x u - KU)Mp e ) + 0(e). 
Let g e — * g° as e — + 0, then it is clear that 

Q(g°)=P- V X /9°M + lu ■ Kp°M + ( t p°u(x) ■ V P M + p a (( r n x (V x u ■ n) - KU) ■ V W M. 

Again note that Q(/) is a linear operator, we may assume the ansatz for g° as 
follows. 

g° = a ■ V x p° + b ■ Kp° + ( t p°c ■ u(x) + p°(C r n x V x u ■ n - KU) ■ d. 

In comparison with the expression of Q(g°) above, we see that it is sufficient to find 
a,b,c,d such as 

Q(a) = pM, 
Q{b) = uM, 
Q(c) = VpM, 
Q(d) = V U M. 

In order to uniquely determine {a, 6, c, d} we state the following 
Lemma 3.5. Let g G L 2 (dpd n uj) such that J gdpd n oj = 0. The problem 

QW=g, (3.2i) 

has a unique weak solution in the space iJ 1 (c?pd n w). 

Proof. For each fixed n G § 2 , we apply the Lax-Milgram theorem to the following 
variational formulation of (|3.21l) : 

k B T J [&MVpV • Vpc/> + £ r MV W V • V w 0] dpd n u = J gcj)dpd n uj (3.22) 

for all <p G H l (dpd n uj). The function M is bounded from above and below on 
IRp x T n (S 2 ), so the bilinear form at the left-hand side is continuous and coercive on 
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H 1 (dpd n u>). The right-hand side is a continuous linear form on H 1 (dpd n u>) due to 
the zero average of g over M.^ x T n (S 2 ) . □ 

Lemma 13.51 ensures that there exists unique a, b, c, d in the space {<p\ J <pdpd n Lj = 
0} since 



J {p M,.M,V t M,V„M }d ^ = 0. 



Note also that from 



V P M = —r-=pM, V^M = —r^—wM, 
k B T k B T 

it follows 

a , b 



k B l k B l 



It remains to determine only a and b. Let a = aMp we have 

Q(a) = -a( t pM = pM, 



which gives a — — . leading to 



Ct 

Similarly let b = (3Moj, then 

Q(b) = ( r k B TV u ■ (PM) 
which together with Q(b) = Mu> gives f} = — 1/Cr, thus 



1 pM. (3.24) 



Cr 

It is straightforward to verify that 



b = -— ujM. (3.25) 



{6, d} ■ pdpd n uj = 0, / {a, c} • ujdpd n uj = 0. 
These enable us to determine the asymptotic limits of the fluxes. 
J i = j pg°dpd n uj 

= J (a ■ W xP °)p + J(b- Rp°)p + CtP° J (c • u(x))p + p° j d- (C r n ■ V x u - KU)p 
= a - \V x p I pdpd n uj, 

J® = J ujg°dpd n uj 

(a ■ V xP °)lu + ( b ■ R P °uj + ( t p° I (c • u(x))uj + p Id - (( r n ■ W x u — 1ZU)oj 



( P° \ 

b ■ 7Zp° — - — —(Crn ■ \7 x u — 1ZU) I ujdpd n u>. 
V k B T J 

In order to further simplify the above fluxes, we state the following lemma: 
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Lemma 3.6. Consider the space pel 3 and (n,u>) G TS> 2 . For any vector A we 
have 



J M(p ■ A)pdpd n uj = kgTA J Mdpd n u, 
J M(u) ■ A)ujdpd n uj = k B T(Id -n® n)A / 

Apply this lemma to the above expressions to obtain 

J\ = lp°u(x) - ^V x p°] J Mdpdnu,, 



(3.26) 
(3.27) 



J° = (Id - n® n) 



-^Tlp + ^-(( r nx V x u ■ n - KU) 



(3.28) 
Mdpd n uj. (3.29) 



A simple calculation shows for arbitrary A that 

(Id — n ® n)(n x A) = n x A. 

Therefore 



J? 



Cr 



n x W x u ■ n TZU 

Cr 



Mdpd n uj. 



Now the limiting equation (|3.20|) divided by C — J Mdpd n w reduces to 

knT . n kuT 



d t p° + V x ■ (u(x)p°) + K-((nxV x u-n~ ^)p°) = ^-A x p° + ^-K ■ Kp° 

Qr 



Cr 



Kip 



k B T' 



■nu 



Regrouping with p° replaced by ip leads to 

d t ip + V x ■ (u(x)tp) + TZ-(nx V x u ■ ip) = D t A x ip + D r TZ 
where 

D t = ^, D r = k -f. 

Ct Cr 

This is exactly the equation derived for rod-like polymers with no incrtial force. 



Proof of Lemma I3.6I The integration is over to from a tangent bundle with 
(n,ui) G T8 2 . We fix n G § 2 arbitrarily, it suffices to evaluate 



2 / a (A ■ u)ujd n oj = ^(Id - n ® n)A I e^^d, 



U). 



Let K be a rotational operator in M 3 such that 

n = Ke 3 , K G 50(3), 
and use transform to = Kui', satisfying \uj\ 2 = \oj'\ 2 . 
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Let A = KB, the transformation gives 

e- u>2/a {A-io)u J d n u J = { e~ [uJ ' )2/a {A- Kw')Ku'duj' 



n-uj— 



J uj' 3 — 




*/. 


e 


^3 


=0 




&i 
















~K(B - e 3 b 3 ) I e~^/ a d n 




e -«^) 2 + (^) 2 )/a duJ l idul l 2 



uJ 



-{Id~n®n)A \ e~" /a d n ui. 



n-uj—0 



Here we have used the fact that 



e 3 & 3 = K 1 nb 3 = nej K T A = K 1 n® nA. 

4. Concluding remarks 

In this work we have derived novel kinetic equations for Dumbbell-like polymers 
as well as rod-like polymers. Inertial forces are taken care of by an augmented 
environment in an extended configuration space. In the case of rod-like polymers, 
the augmented space for orientation is just a tangent bundle of the usual sphere. 
We have also shown that the formal limit of the augmented equation recovers the 
usual inertia-free kinetic models explored in literature. 
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